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ABSTRACT 

We present analytic models for the formation and evolution of tidal tails and related 
structures following single or multiple impulsive disturbances in galaxy collisions. Since 
the epicyclic approximation is not valid for large radial excursions, we use orbital 
equations of the form we call p-ellipses (a class of precessing ellipses). These have 
been shown to provide accurate representations of orbits in logarithmic and power- 
law halo potentials. 

In the simplest case of an impulsive collision yielding a purely tidal disturbance 
the resulting tidal tails have simple structure. Scalings for their maximum lengths and 
other characteristics as nonlinear functions of the tidal amplitude and the exponent 
of the power-law potentials are described. The analytic model shows that azimuthal 
caustics (orbit crossing zones of high density also seen in numerical models) are pro- 
duced generically in these tails at a fixed azimuth relative to the point of closest 
approach. Long tails, with high order caustics at their base, and ocular waveforms are 
also produced at larger amplitudes. 

The analysis is then extended to nonlinear disturbances and multiple encounters, 
which break the symmetries of purely tidal perturbations. The p-ellipse orbital so- 
lutions are similar to those in the linear tidal case. However, as the strength of the 
nonlinear terms is varied the structure of the resulting forms varies from symmet- 
ric tails to one-armed plumes. Cases with two or more impulse disturbances are also 
considered as the simplest analytic models distinguishing between prograde and retro- 
grade encounters. The model shows explicitly how tail growth differs in the two cases. 
In the prograde case a specific mechanism for the formation of tidal dwarf galaxies 
at the end of tails is suggested as a consequence of resonance effects in multiple or 
prolonged encounters. Qualitative comparisons to Arp Atlas systems suggest that the 
limiting analytic cases are realized in real systems. For example, we identify a few 
Arp systems which have multiple tidal strands meeting near the base of long tails. 
These may be swallowtail caustics, where dissipative gas streams are converging and 
triggering star formation. UV and optical images reveal luminous knots of young stars 
at these 'hinge clump' locations. 
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1 INTRODUCTION 

Spiral waves, and their longer extensions which we now rec- 
ognize as tidal tails, have fascinated since the studies of 
Holmbergl fl940| |T9ll) , |Zwicky| ( |T959| ) and others. Holm- 
berg and Zwicky speculated that tails, and related struc- 
tures, were produced by tidal forces in galaxy collisions. 
This was well confirmed by the numerical case studies de- 



scribed by Toomre & Toomre ( 1972 with references to other 



modeling efforts). These results have in turn been amply 
confirmed by subsequent generations of increasingly com- 
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prehensive numerical simulations, including detailed studies 
of specific systems (see e.g., the recent papers of Struck, et 
aT1[iu05l [Smith, et al.|[2008| |Chien fc Barnesl|2010[ |Barnes 



2011 



and references therein), and studies of detailed fea- 
tures that can develop within the tail, like tidal dwarf galax- 



ies (e.g. 


, Barnes & Hernquist|1992b 


Elmegreen, Kaufman & 


Thomasson 1993, Due, Bournaud & Massct 


2004 


Bournaud 


& Due 


2006 


Bournaud, Due & Emsellem 


2008 


. The now 



extensive literature of numerical models provides a broad 
understanding of the basic processes of tail formation and 
evolution. However, as usual with simulations, it is hard to 
grasp all the systematics of a complex phenomenon like tidal 
tails even with the results of many of specific calculations 
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available. Analytic models or scalings are often very helpful 
for understanding such trends. 

A number of papers in the literature have used ana- 
lytic or semi-analytic models to study waves produced in 
galaxy collisions. One of the first of these was Toomre's (see 
Toomre|1978} also |Struck-Marcell fe Lotan|1990| illustration 
of the radial orbital motions that make the circular waves of 
a colliding ring galaxy. Later came a number of analytical at- 
tempts to map out less symmetric evolving wave structures 
through disturbed discs using the Impulse Approximation 
for the perturbation, and the epicyclic approximation for 



initially circular orbits following the perturbation (Struck 
Marcell|1990| |Donner, Engstrom fc SundeliusffTWTj [Gerber 



fc Lamb|1994 1. In these studies the orbits are approximated 



analytically, but the waves were studied by numerical calcu- 
lations of the post-collision orbital motions of sample par- 
ticles in a model disc. The degree to which such models 
can match more sophisticated, self-consistent simulations is 
generally impressive at early stages. However, these initial 
collisions are often very impulsive, and the full simulations 
are required to model detailed dynamics of the final merg- 
ers of galaxies (e.g., reviews of Barnes & Hernquist 1992a 



Struck|1999 1. This limitation of the semi-analytic models is 



not surprising. More vexing was the fact that the evolving 
wave structures are complex, and seem to depend sensitively 
on values of collision parameters. Thus, it was not immedi- 
ately obvious that the semi-analytic models could provide 
any more predictive power than case-by-case studies with 
simple simulations, e.g., restricted three-body calculations. 

However, there is perhaps more hope for analytic mod- 
els of the simplest and most symmetric galaxy collisions, as 
in the case of symmetric ring galaxies. In a series of papers, 
one of us has derived a complete analytic model of stellar 
ring waves, subject to the Impulse and epicyclic approxi- 
mations, and the further approximation of rigid halo poten- 



tials ( jAppleton fc Struck-Marcell||1996| |Struck| |2010). The 

key to this work was the realization that the orbit-crossing 
zones that define the waves can be described with the math- 



ematical machinery of caustic or catastrophe theory ( Arnold 
19861). 



More recently, we have been able to extend this analytic 
theory to cases of spiral waves produced in impulsive flyby 
encounters ( jStruck, Dobbs fc~ Hwang 2011). This theory ap- 
plies in cases where the galaxy discs have a relatively high 
value of the Toomre Q parameter, so that the induced spirals 
are not strongly self-gravitating, and the epicyclic approx- 
imation can be used to describe the post-encounter orbits. 
The obvious difference from the ring galaxies is that the 
magnitude and direction of the velocity impulse, and thus, 
the parameters describing the epicycles, vary with azimuth. 
The radially and azimuthally dependent velocity impulses 
can be simply approximated by classical tidal formulae, or 
more accurately in specific cases by formulae like those de- 



rived, for example, by D'Onghia, et al. (20101. 



The first difficulty that confronts anyone wanting to ex- 
tend such analytical models to tidal tails is their length. 
Tidal tails can stretch out to distances beyond twice the orig- 
inal orbital radius of the stars or gas clouds within them. Us- 
ing a single epicycle to approximate disturbed orbits works 
best when the disturbance is relatively small, and the radius 
of the epicycle is a fraction of the original, or the guiding cen- 
tre, radius. Thus, the epicycle approximation applied to tail 



particles is an absurdity. In such encounters, stars in two disc 
quadrants have their angular momentum increased, and fly 
out in bridges or tails. Those in the other two quadrants lose 
comparable amounts of angular momentum, and in simula- 
tions these stars generically pile up in the central disc. The 
epicyclic approximation cannot represent this phenomenon 
well either. 

Another problem in analytically modeling tails is that 
the conventional wisdom has it that long tails cannot be 
produced in impulsive encounters; more torque is needed 
to pull tail particles far from the disc centre. Specifically, 
the Impulse Approximation does not distinguish between 
prograde and retrograde encounters, and D'O nghia, et al7| 
( |2010[ ) argue that it is nearly resonant encounters in the for- 
mer case that produce long tails. These authors compute the 
perturbation to circular stellar orbits due to the finite-time 
pull of a companion on a fixed orbit. They suggest that one 
can assume that these (for some stars semi-resonant) per- 
turbations are produced instantly in a kind of semi-impulse 
approximation, which would distinguish between prograde 
and retrograde encounters. This work has similarities to the 
Linear Response Theory approach of |Nelson fe Tremaine| 
fl999| ). Both provide means of improving the classical Im- 
pulse Approximation. In the present work we will use the 
classical approximation, with some extensions. We will also 
probe the limits of the generalization that tails cannot be 
produced in impulsive encounters. 



But how should we produce such analytic models with- 
out the simple and versatile epicyclic orbit approximation? 
We are forced to look at other types of orbits, and a variety 
of other forms have been suggested, even in recent years. All 
of them ultimately have the approximate form of precessing 
ellipses. However, |Struck| ( |2006| henc eforth S06), |Lynden-| 
Bell fe Jin| (|2008[ henceforth LBJ) and|Lynden-Bell| (|2010a|) 



rediscovered a particular class, whose members are described 
by simple formulae, and which provide good fits to orbits in 
power-law and logarithmic potentials, ranging from nearly 
circular to nearly radial. The fits are very good for near 
circular orbits. They allow errors at second-order in the ec- 
centricity, but since they accurately approximate the preces- 
sion period, they do not drift in phase relative to the true 
orbit. (Better, but more complex approximations were also 
explored in S06.) The symmetry properties of these curves, 
which we call p-ellipses (for precessing ellipses in power-law 
potentials), have been explored by LBJ. The Impulse Ap- 
proximation, and orbits described by these curves provide 
the basis for constructing analytic models of tails, as we will 
shortly describe (also see Lynden-Bcll 2010b I. 



Before beginning that discussion we note that we will 
restrict attention to free tails without connection to a com- 
panion. Bridges experience prolonged forcing from the com- 
panion, so we cannot expect the Impulse approximation to 
apply. The models discussed below do not include the effects 
of self-gravity in disc waves, and are two-dimensional. That 
is, the orbital plane of the companion is assumed to be the 
same as the disc plane of the target, and there is no vertical 
perturbation. It would be straightforward to generalize the 
models to include moderate vertical effects, but we do not 
do so in this paper. 
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2 ANALYTIC MODEL EQUATIONS AND 
THEIR PROPERTIES 

2.1 Post-impulse p-ellipse orbits 

To begin the analytic study we assume, for simplicity, that 
the disc orbits are circular before the encounter. The polar 
equation for p-ellipse orbits can be written, 



£ = [1 + ecos (mfl)]5 +i . 



(1) 



where the orbit of an individual star is specified by the semi- 
latus rectum p, which sets the scale of the orbit, the eccen- 
tricity e, and an initial value of the phase 9 a . The parameter 
S is determined by the form of the gravitational potential, 



e.g.. 



for the logarithmic potential, and k for the Ke- 



pler potential. The factor m is a simple function of S, see 
S06 for details. (In S06 this factor was written as (1 — 6), 
but here we use the simpler notation of LBJ.) 

Assuming that the form of the potential is fixed, after 
the impulse the problem of determining the p-ellipse orbit 
of a star is equivalent to finding the values of p, e, and 8 
(which determines the initial position of the star on the p- 
ellipse). A star's position in the disc is specified by radius 
r, and azimuth, <f). We denote the coordinates of a star at 
the time of impact as r — q, and <f> = cj> . In the Impulse 
Approximation these are unchanged immediately after the 
impact; only radial and azimuthal velocities are changed. We 
choose the coordinate system such that the impact occurs 
at azimuth <j> — 0. Equation (1) shows that when the p- 
ellipse azimuth — 0, r takes its minimum value, which 
does not generally occur at <j> = 0. Thus, the two azimuthal 
coordinates are offset, and in general o 7^ O for any star. 
Note that we use the symbol 9 for the p-ellipse azimuth in 
this work, while <j> was used in S06. 

Let us return to the issue of finding relations that can 
be solved for the values of p, e, and 9 for a stellar orbit. 
The first is given by equation (1), written in terms of initial 
values, 



[1 + e cos (m0 o 



li+<5 



(2) 



To obtain a second relation we differentiate equation (1) to 
get a radial velocity equation. This can be written as, 



pr 



(|+5) mesin {m0)0 
[l + ecos(m6>)]5~' 5 ' 



(3) 



where r and 8 are the time derivatives of the coordinates, 
i.e., the radial and angular velocities. The initial values of 
these velocities are determined by the Impulse Approxima- 
tion. Equation (1) can be used to eliminate p/r, and equa- 
tion (3) can be solved for the ratio of velocities, and evalu- 
ated at the time of impact as a function of q and Q a . I.e., 



v r (| + 5) mesin (m£ 
vq [1 + e cos (m0)] 



(4) 



where v r — r, and vg = qO. Note that due to fortuitous 
cancellations there are no exponents with factors of 8 in 
equation (4). 



Equation (4) is the second relation between the p-ellipse 
parameters. A third relation can be obtained from the con- 
servation of angular momentum, and these equations can 
then be solved for the three orbital parameters. It turns 
out, however, the symmetries of simple tidal perturbations 
allow us to avoid this latter calculation. 



2.2 Simple tidal perturbations 

As noted above, the encounter perturbation can be approx- 
imated at various levels of accuracy. iD'Onghia, et al. ( 2010 1 



give some of the most sophisticated analytic formulae pub- 
lished to date, and they are quite complex. |Gerber fc Lamb] 
( 1994 1 computed the nonlinear impulse from a Plummer po- 



tential perturber on a parabolic orbit, which is less complex, 
though still not simple. When combined with the equations 
for the p-ellipse orbits, complex expressions for the perturba- 
tion velocities can render further calculations very difficult. 
To begin to explore what can be learned from analytic mod- 
els, we begin with the lowest order impulse formulae derived 
from equation (7) of Gerber and Lamb. We call this the Pure 
Tidal case. These impulses can be written, 



Av r = A 



V^o cos (2</> ) , 



v^o sin (2<j)o) 



with, 



A = 



2GM C 

1 rain 



(5) 



(6) 



(7) 



In these expressions r m in is the distance of closest approach 
between the galaxy centres, V is the relative velocity at 
closest approach, M c is the total mass of the companion. 
Equations (5) and (6) are the approximation of Gerber and 
Lamb's equation to first order in the quantity q/r m in, which 
are nonlinear in this ratio (see Section 5). Note that these 
expressions differ from textbook forms for the tidal acceler- 
ations in the Earth-Moon system and other published forms 
(e.g., the factor of 2 in the argument of sine and cosine). 
These differences trace to different assumptions about the 
acceleration of the centre of the disc galaxy. 



2.3 Orbital parameters 

Equations (5) - (7) can be combined to obtain an expression 
for the post-impulse velocity ratio, analogous to equation 
(4), 



V r 



Av r 



A' cos (20 o ) 



v<{, + IS.V4, 1 — A' sin (2<p ) ' 
In this equation we use the radius-dependent amplitude 



(8) 



A' — A 



(9) 



The form of equations (4) and (8) is very similar. As dis- 
cussed in Sec. 2.1, we expect to fit a post-collision orbit 
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with a p-ellipse. We can do this by setting equations (4) and 
(8) equal to each other, 



(| + 8) me sin (m8 ) 
[1 + e cos (m0 o )] 



A' cos (2<f> ) 
- A' sin (20 o ) 



(10) 



and requiring that the numerators and denominators (radial 
and azimuthal velocities) on each side are equal. 
The solution is obtained by taking, 



m0 o 



2<)> + - + ip, 



(11) 



where the 7r/2 term turns sine into cosine and vice versa, 
and the angle ip gives the remaining phase difference. With 
this relation between the angles, we equate the numerators 
and denominators of equation (10), and solve for A' and xp. 
We get, 



A' = 



sf{\ + S) 2 m 2 sin 2 (2<t> ) + cos 2 (20 o ) ' 



(12) 



and, 



tan (20 o + i/;) = ( i + 5 J m tan (2<j> ) = 

I + Sj \/2(l-5) tan(20 o ). (13) 



where the last equality in equation (13) makes use of the first 
order approximation from S06 for m(8). Equation (12) gives 
A' ~ e for most values of O . Thus, the eccentricity of the p- 
ellipse depends weakly on the initial azimuth. The equations 
also show that perturbation amplitudes of order unity are 
required to get high eccentricities, and nearly radial orbits. 

The solution to equation (10) given by equations (11)- 
(13) is very simple. Note that the use of equation (11) in 
equation (1) also yields an expression for the semi-latus rec- 
tum of the p-ellipse, 



1 + e 



cos \ 2<f)o + — + ip 



l>\ = 1- 



' JT\ + S) 2 w? sin 2 (2^0 + cos 2 (20 o ) \ 



(§ + *) 



(14) 



This parameter depends on both the star's initial radius and 
azimuth. 



2.4 Accuracy of the approximation 

The obvious question about the simple approximate orbits of 
the previous subsection is how accurate are they, even within 
the assumed fixed potential of the parent galaxy? The imme- 
diate answer is that p-ellipse orbits in a power-law potential 
have been shown to be quite accurate (S06, LBJ, |Lynden^| 
Bell||2010a l. Though initiated by an impulsive disturbance, 
the orbits described above will inherit those general charac- 
teristics. For example, one of the most remarkable properties 
of p-ellipse orbit approximations is that in many cases the 




Azimuthal Advance 



Figure 1. Radius versus azimuthal increase from the initial value 
(as a proxy for time) is shown for sample numerically integrated 
(solid curves) and analytically approximated (dashed curves) stel- 
lar trajectories. The initial radius in all cases is q = 1.0, but suc- 
cessive pairs of curves are offset by one unit in radius for graphical 
clarity. Successive pairs of curves are distinguished by different 
initial azimuths relative to the line connecting the disc centre to 
the point of closest approach (i.e., the impulse direction). From 
bottom to top the initial azimuths are: 0.3, 0.6, 0.9, 1.2, 1.5, 1.8, 
2.1, 2.4, 2.7, 3.0, 3.3 radians. The form of the potential is logarith- 
mic (<5 = 0) and the amplitude of the tidal impulse is A' = 0.50 
(at q = 1.0) in all cases. In the analytic approximation the eccen- 
tricities were computed using equation (12). 



fit to the true orbit does not deteriorate with time. That 
property also continues to hold. On the other hand, while 
the above p-ellipse approximation is quite accurate for low 
to moderate eccentricity orbits, it is significantly less accu- 
rate at higher eccentricities. Tidal tails are made of quite 
eccentric orbits, so we will consider this issue more carefully 
below. 

We begin that discussion with Figure 1, which com- 
pares 11 numerically computed orbits, each at a different 
initial azimuth relative to the impulse impact azimuth, and 
approximate orbits given by the equations of the previous 
subsection in the logarithmic potential (5 = 0). The initial 
radius for all the orbits shown is the outer disc radius, which 
is scaled to 1.0. In the graph, successive orbit pairs are offset 
by one unit for clarity. The impulse amplitude in the case 
shown in Figure 1 has a relatively high value of A' = 0.5. 
Nonetheless, the numerical orbits (solid curves) are quite 
well fit by the analytic approximation (dashed curves). The 
maximum deviations in radius between the two curves are 
about 50% (of the larger radius), and generally much less. 

Within this figure it can be seen that the fit is better for 
some trajectories than for others. The fit correlates directly 
with eccentricity, which correlates with amplitude according 
to equation (12), but also with initial azimuth via the tidal 
perturbation. Higher values of the eccentricity in the post- 
impulse orbit yield worse fits. The case shown in Figure 1 is 
in fact a borderline case. The fit is generally much better for 
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Figure 2. Same as Figure 1, but with an impulsive tidal ampli- 
tude of A' = 0.80 (at q = 1.0). 

low amplitude cases, and gets worse rapidly as the amplitude 
increases above 0.5. This can be seen in Figure 2, which 
shows a higher amplitude case. This is easy to understand. 
The analytic solution solves the equation of motion to order 
e. Beyond A' , e ~ 0.5, the effect of the nonlinear terms grows 
rapidly. 

The fact that a simple p-ellipse approximation works 
up to amplitudes as high as 0.5 is impressive. At somewhat 
higher amplitudes, the eccentricities of specific orbits can 
exceed values of 0.9, and the difference between numerical 
and analytic curves can exceed factors of a few. It is true 
that only a small fraction of all orbits are affected, and they 
could be neglected, but this fraction grows with amplitude. 
Thus, the approximation becomes increasingly less useful. 

Fortunately, a little further analysis allows us to largely 
overcome this limitation. The inaccurate analytic eccentrici- 
ties occur when the cosine term in equation (12) is large. An 
ad hoc, but modest change of that equation, to the form, 



J (§ + S) 2 m 2 sin 2 (20 o ) + (jhj) cos2 ( 2< M 

proves very helpful. The added factor of 1/(1 + A') largely 
cancels the effects of the nonlinear terms, and the modified 
analytic orbits fit the numerical ones quite well. This is true 
throughout the range 0.5 < A' < 0.9, while the effect of 
this extra factor is negligible at small amplitudes. With this 
modification we have a qualitatively good orbit approxima- 
tion over most of the interesting range of amplitude. 

The examples shown in Figures 1 and 2 are for the 8 = 
case. The improved approximation of equation (15) does not 
work as well for larger (or smaller, negative) values of 8. In 
the Keplerian case (8 = 0.5), it starts to break down for 
amplitudes greater than about 1/3. Specifically, the maxi- 
mum radial excursions of some of the numerical orbits grow 
very large as the amplitude increases. This is not a surprise. 



Figure 3. Same as Figure 2, except that in the analytic ap- 
proximation the eccentricities were computed using equation (16). 
Comparison with Figure 2 shows the improvement, see text. 

An azimuthal pertubation amplitude of 2 1 / 2 — 1 boosts the 
stellar velocity to the classical escape velocity. There is no 
such problem in the logarithmic potential, but there will be 
in all the potentials with values of 8 greater than zero. 

Equation (15) is designed to tame the excessive excur- 
sions in the analytic orbits at high amplitude in the 8 = 
case, and a similar modification of the sine term in that equa- 
tion can boost the analytic solution at the correct phases to 
match the larger numerical excursions in the Keplerian and 
intermediate cases. For example, the form, 

A' = 

(j+j)me 

^b 2 (| +5) 2 m 2 sin 2 (20 o ) + b! cos 2 (20 o ) 

with, bl =( r ^ 7 ), b2 = (l + A')°'™ 2 , (16) 

yields fits better than or comparable to those of Figures 1 
and 2 for all values of 8 between and 0.5, and for all, except 
the highest values of amplitude. This is shown in Figure 3, 
which is the same as Figure 2 except that the eccentricities 
are derived from equation (16) rather than equation (12). 
The effect is even more pronounced for cases when 8 > 0. 
The new coefficient of the sine term was derived by experi- 
mentation. It has not been proven to be an optimal fit, but 
it is a good one. The 0.5 term in the exponent of 62 off- 
sets much of the difference between analytic and numerical 
curves at 8 = (see Figs. 1 and 2). A nonlinear dependence 
on 8 is needed in the exponent to fit the range of potentials, 
and the large coefficient is needed to offset the low numerical 
value of 8 2 . 

It is worth recalling that the refinements of the simple 
orbit approximation of the previous subsection embodied in 
equations (15) and (16) only change the value of the eccen- 
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tricity. They still yield a p-ellipse fit to a given orbit, but one 
that is better than the fit achieved by the velocity matching 
of equation (10) at certain azimuths. That velocity match- 
ing gives a p-ellipse solution that is guaranteed to fit the 
orbit to first-order in e. In sum, we have shown that there is 
also a p-ellipse, with a different eccentricity, that fits some 
perturbed orbits better. In the present context this may ap- 
pear fortuitous, but in fact it is a corollary of the result that 
any orbit in these power-law potentials can be well fit by a 
p-ellipse (S06, LBJ). We also note that D. Lynden-Bell has 
explored extended versions of the p-ellipse functions that 
provide still better approximations at high eccentricity (pri- 
vate communication). 

2.5 Special case: the backbone of the tail 

Tails are produced by the strongest azimuthal impulses, and 
equation (6) implies that these will occur at azimuths of 
(f> = — 7r/4 and 37r/4, relative to the line connecting the 
disc centre to the point of closest approach. The second of 
these azimuths is the one associated with the counter-tail. 
The radial velocity impulse is zero at these azimuths. 

Of special interest is the structure of the p-ellipse orbits 
of stars initially located along this azimuth, which can be 
regarded as the backbone of the tail. Since the radial impulse 
is zero at this azimuth, equation (4) yields sin(m(9 ) = 0, 
m8 = 0, or 2-n, and ijj = according to equation (11). 
Substituting this value into equation (14) we get, 



P 



= (1 + A') 



(*+*) 



(17) 



Thus, in going from values of S = (logarithmic) to 0.5 (Ke- 
plerian) , the orbit scale p goes from a weak square root to a 
linear dependence on A' . The difference is even greater if we 
use equation (16), with its extra terms, for the eccentricity 
in equation (14). 

Using equations (1), (14) and (16), we obtain the follow- 
ing expression for the maximum radial excursion (apoapse) 
of a backbone star, 



1 + A' 
1-e 



1 + A' 



A' (1 + A')°- 25+1 



(*+*) 



(18) 



This equation has several interesting properties. The first 
is that the denominator goes to zero at a critical value of 
the amplitude A' ^ 1, so the expression allows for very long 
tails. If the colliding galaxies are bound, it can be shown that 
the amplitude is equal to a factor of order unity times the 
ratio qM c /(r m i„M(q)), where M c is the mass of the com- 
panion, and M(q) the mass of the parent galaxy within the 
star's orbit. If the two masses are comparable and the dis- 
tance of closest approach not much greater than q, then the 
critical amplitude can be exceeded. E.g., a close encounter 
with a more massive (and compact) companion will produce 
strong perturbations, and according to these results could 
send backbone stars out to large distances. When S > stars 
can escape the potential as noted above. (These conditions 
are apparently satisfied in Arp 105 and 174, for example, not 



to mention the globular cluster Palomar 10.) This situation 
would probably entail the near destruction of the smaller 
galaxy, but it does suggest that in some cases, impulsive 
encounters can produce long tails. 

In the case of weaker disturbances, equation (18) sug- 
gests that it would be hard to launch a backbone star out 
to more than a few times its initial radius. For example, 
with 5 = (logarithmic potential), and A' = 0.3, then 
r ap /q ~ 1-4. 

The dependence of equation (18) on the form of the po- 
tential is also interesting. If we take 8 = 0.5 and A' = 0.3, 
then r ap /q ~ 4.4. Other factors being equal that implies 
much greater stretching of tails in Keplerian potentials than 
in flat rotation curve potentials (see Sec. 3). As the ampli- 
tude is increased to A' = 0.4, the tail length ratio increases 
greatly. [Due, Bournaud fc Ma sset (2004!) noted this effect in 
their large grid of tail-making galaxy collisions. They argued 
that magnitude of the stretching was a key factor in deter- 
mining whether a tidal dwarf galaxy could form at the end 
of a tail. If the tail material is stretched too much it is much 
more difficult for local self-gravity to assemble a dwarf. The 
earlier simulations of Mihos, Dubinski & Hernquist ( 1998 1 



also showed that dwarfs are easier to manufacture in more 
extensive halos. The equations above suggest scalings for 
the tail stretching, and the basis of a method for estimating 
other tail properties. 

Another aspect of the tidal stretching results from the 
radius dependence of the amplitude A'. In equation (17) 
this nonlinear dependence results in greater stretching of 
stars with larger initial radii. This global stretching means 
that kinematic pile-ups and radial caustics are not possi- 
ble in impulsive tails, but must be the result of prolonged 
disturbances or local self-gravity. 

2.6 Fleshing out the backbone 

In this section we briefly demonstrate that post-collision or- 
bits in nearby sectors are very similar to those on the back- 
bone. Qualitatively we expect this, given the smoothness 
and continuity of the simple tidal perturbation, but it can 
be demonstrated more quantitatively. We begin with equa- 
tion (12) for the eccentricity. (The results here will also hold 
for the more complex form of equation (16).) For sectors 
near the backbone, we assume that <f> m 3n/4 + e, where 
e << 1. Substituting this into equation (12), solving for e, 
and expanding in powers of e, we obtain, 



e = A' 



2e' 



1 - 



(19) 



Thus, the eccentricity has only a relatively weak, second- 
order dependence on e. 

A similar expansion of equation (15) shows that the 
semi-latus rectum also depends weakly on the azimuthal de- 
viation e, 



.\'>m[ 37r +2e 



(l + A')^ +s 



2c 



1 + A' 



(20) 
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The weak dependence on e in both p-ellipse orbit parameters 
confirms our statement that adjacent orbits are like those in 
the backbone, and that the backbone represents the whole 
tail. 



2.7 A sharp wave in the tail 

Section 2.5 ended with the comment that tidal stretching 
rendered kinematic crowding and caustic wavefronts unlikely 
in the radial direction. This is not necessarily the case in 
the azimuthal direction. Technically, the caustic condition 
requires that the Jacobean matrix (which measures com- 
pression and shear) of the coordinates with respect to their 
initial values go to zero. We assume that the derivative of 
the radius with initial azimuth (dr/d(j) ) is negligible, i.e., 
nearly uniform stretching in azimuths near <j> — 37r/4, in 
accord with the results of the previous section. Then the 
condition for caustic compression reduces to the requirement 
that d<j>/d<t>o = 0. 

An expression for the time dependence of the azimuth 
of a p-ellipse orbit is given by equation (C2) of S06, which 
can be written, 



tan 



1 + e 



tan 



tan 



m0 o 



(21) 



where h is the specific angular momentum of the orbit and 
t is the time since impact. In general the parameters p and 
e will depend on the initial azimuth, (f> . However, in accor- 
dance with the previous section, the azimuthal dependence 
of e and p is weak at azimuths near <f> = 3n/4. Thus, when 
differentiating equation (21) with respect to cj> , we neglect 
the derivative of these parameters. 



d 



tan 



m9 



tan 



m8 



0. 



(22) 



Combining this with the condition for the appearance of a 
(azimuthal) caustic we find, 



= = 



W 



= 2 (mBo 



This is satisfied when, 



= - = e + A<t>. 

m 



(23) 



(24) 



For stars initially on the backbone this reduces to 9 = 7r/v2 
and (j> — 5n/4, when 5 = 0. This value of 6 shows that a star 
encounters the caustic at a fixed phase (outward moving) of 
its p-ellipse orbit. Since the tail is primarily made up of stars 
originating near the backbone the second equality implies 
that the azimuth of the tail caustic is close to (j> — 57r/4. At 
this level of approximation the shape of the tail is close to 
a straight line, while the caustic persists. 

Like the caustic edges of collisional ring waves discussed 
Appleton fc Struck-Marcell| ( |1996} also |Struck|2010[ |, this 



tail caustic is an archetype. It is a generic consequence of 
any flyby disturbance strong enough to produce a clear tail. 



Equation (24) only gives a leading edge caustic. Equa- 
tion (23) can give others (e.g., a trailing caustic), and a less 
approximate treatment of the caustic condition could yield 
still more. The tail caustics generally end in an outer cusp. 
The inner base could either be another cusp, or a more com- 
plex caustic. Numerical images described below show that 
swallowtail caustics (see e.g., Arnold|1986 1 are common. 

The leading edge caustic is commonly seen in self- 
consistent numerical simulations of tails (e.g., Gerber fc 
Lamb 1994). Given that gas will pile up there, it is a likely 
site for star formation. This is especially true at the early 
stages of tail development when strong compression in the 
caustic, and as yet moderate stretching along the tail, may 
give optimal conditions. 



2.8 More caustics and swallowtails 

The backbone sector represents the basic tail dynamics, but 
the sectors 90° away at initial azimuths cj> = tt/4 and 57r/4 
tell a different story. According to equation (11), at these 
azimuths m0 o = ir and ip = 0. Equation (14) yields, 



P 
<1 



(1 



,i+s 



The analog of equation (18) in this case is, 



q 



1 



I -A' 
1-A' 



(*+') 



1. 



(25) 



(26) 



Thus, the orbit size is decreased by the impulse along these 
sectors, and the initial radius is the apoapse. The radial ve- 
locity impulse is zero along these sectors and the azimuthal 
impulse is opposite the star's rotation. Stars along these and 
nearby sectors will fall inward and pile up at periapse, form- 
ing caustics in the inner disc. Alternately, at large radii, 
where the perturbation is strongest, they can contribute to 
the intersecting stellar streams that make swallowtail caus- 
tics at the base of the tidal tail. 

The sectors around initial azimuths of (f> = and 7r help 
complete the story. They are pulled outward, away from the 
nucleus by the tidal perturbation. That is, the disturbance 
provides an outward radial velocity, and small or zero az- 
imuthal velocity perturbation. This perturbation is like that 
in the disc of a colliding ring galaxy, except 180° out of phase 
(pulled outwards instead of inwards). Stars in these sectors 
will first move outward, and then inward. In the outer disc 
they can contribute another star stream to the swallowtail 
caustics at the base of the tidal tails. Swallowtails have five 
intersecting streams, in contrast to the three in ring wave 
caustics. 

Stars in sectors around initial azimuths of <f> = n/2 and 
37r/2 will be pulled inward by the tidal perturbation even 
more like ring galaxies. However, this pull is relatively weak, 
and at best only generates a weak wave in the outer disc. 



2.9 Ocular waves 

It is worth briefly examining the approximate radial caus- 
tic condition before leaving the topic of caustic waves. This 
condition is dr/dq = 0, i.e., the complete compression of 
an initially finite range of radii. Using equation (2) for p 
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in equation (1), differentiating with respect to q, and then 
simplifying the result we obtain, 



dr 
dq 

1 + 



1 + e cos (m8 ) 



1 + e cos (m6 

de ~ 



+ 5) q 



dq 



cos (m9 ) — cos (m8) 



(1 + e cos (m6)) (1 + e cos (m0 o )) 



(27) 

As long as e < 1, the caustic condition can only be satisfied 
when the second term on the right-hand-side (in the outer 
square brackets) is zero. Thus, the approximate caustic con- 
dition is, 



de 
^ dq 



i (m0 o ) — cos (m0) 



(1 + ecos(m6»)) (1 + ecos(m6» ))J ' ^ 

Note that with e given by equation (12) (or equation (16)), 
with substitutions from equations (7) and (9), the value of 
the factor q(de/dq) will be within a factor of a few of e, 
for the power-law potentials considered here. Given that the 
factors in front of the term in square brackets are less than 
or about unity, we want to minimize the negative value of 
the term within the square brackets to satisfy the equation 
(or make the right-hand-side of the equation negative, in- 
dicating an orbit between two caustics). Evidently, this is 
achieved near m0 o = ir and m9 = 2tt. Equation (12) then 
implies that <f> — 7r/4 or 5n/4 (perpendicular to the back- 
bone). Then, if> is given by, A<f) = A9 (equation (24).) Thus, 
with the values of the other variables just given we find that 
= tt/4 + ""(V 2 + S)m 2 or <f> = 5tt/4 + tt(1/2 + S)m 2 . In 
the logarithmic (Keplerian) potential with 5 = (1/2), then 
4> — 37r/4 or 77r/4 (in both cases) to satisfy the caustic con- 
dition. 

As noted in the previous subsection, stars near these 
initial azimuths lose angular momentum from the impulse 
and fall inward. Equation (28) suggests that their pileup at 
periapse generates a radial caustic if e and its gradient are 
sufficiently large. I.e., when the impulse is sufficiently strong. 
These caustics appear at azimuths orthogonal to the angular 
caustic in the tail. These facts suggest the identification of 
the radial caustics as ocular waves of the type described by 
Elmegreen, et al. ( 1991[ ). The present analysis agrees with 



numerical results in those references in several respects, in- 
cluding the association of ocular waves with nascent tails, 
the orthogonal orientation of these two features, and the ex- 
istence of a threshold perturbation amplitude to generate 
the oculars. Though again very approximate, the present 
analysis suggests that the ocular waves are another caustic 
archetype. 

In the previous subsection we also noted the presence of 
a radial oscillation component if the orbits of stars at initial 
azimuths of (j> — (or n) and <j> — tv/2 (or Sir/ 2). In these 
cases, equation (11) yields values that make cos(m6l ) = 0. 
Then, inspection shows that while the last term in equation 
(28) can easily be negative, it is unlikely to reach a value 
of one and satisfy the caustic condition. Undoubtedly, there 
will be radial compressions of the stars near these initial az- 
imuths, but strong radial caustics are unlikely. We also note 
that slowly developing spiral caustics of the type considered 




Figure 4. Six snapshots of the evolution of an analytic model disc 
in a logarithmic (5 = 0) potential after a tidal impulse at time 
t = 0. The amplitude of the impulse is A' = 0.80 (at q = 1.0), 
the initial disc radius in dimensionless units). The orbital period 
at the initial disc edge is 2tt dimensionless time units. Note the 
different scales of the different panels. See text for details. 



Struck, Dobbs & Hwang (2011 1 are not captured by the 



rough analysis here, though they are related. 



3 GRAPHICAL REALIZATIONS OF THE 
ANALYTIC RESULTS 

3.1 Illustrating the analytical discs and tails 

The simple analytic model of the previous section yields 
many powerful insights on disc and tail structure. However, 
the full tidal structure of this model can only be revealed by 
looking at a large sample of disc orbits described by the an- 
alytic equations. This is best done graphically. Some results 
are shown in Figures 4-8. 

Each of these figures shows the positions of 38,100 stars 
initially distributed uniformly in the disc, and assumed to 
have circular velocities in centrifugal equilibrium with a 
power-law gravitational potential. Subsequently, the orbits 
of these stars were computed according to the analytic solu- 
tions of the previous section following an impulsive distur- 
bance of specified magnitude (and with equation (16) used 
for the eccentricity). In addition the azimuths of the stars as 
a function of time after the impulse on p-ellipse orbits are 
given by a nonlinear, multi- valued function of time (equation 
(C2) of S06). Rather than find the appropriate solutions to 
this equation for every time, we numerically integrated the 
equation for the azimuthal velocity for each star to each 
desired time value. 

In this section we focus on the results for the loga- 
rithmic and Keplerian potentials, which bound the range of 
physically interesting flat-to- declining rotation curve slopes. 
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Figure 5. Same as Figure 4, but with a tidal impulse amplitude 
of A' = 0.50 (at q = 1.0). 
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Figure 6. Same as Figure 4, but with a tidal impulse amplitude 
of A' = 0.30 (at q = 1.0). At this amplitude, or lower, the tails 
appear more like internal disc spirals. 



Specifically, Figures 4-6 each show six snapshots of the post- 
collision disc in the logarithmic potential, with different per- 
turbation amplitudes in each case. 

In these models the units are chosen such that the radius 
of the disc is 1.0, and GM(R), the gravitational constant 
times the mass contained within the disc, is also set to 1.0. 
Then, the circular velocity at the outer radius is 1.0, and 
the orbital period is 2n. In this paper we use only these 
dimensionless units, but a consistent set of scaling factors 
is, for example, R — 10 kpc, v(R) — 250 km s _1 , M(R) — 
1.5 x 10 11 M , and P(R) = 250 Myr. 

Figure 4 shows the highest amplitude case, with A' — 
0.8 for particles at the outer edge of the disc. The first 4 
panels of this figure show the steady growth of the tail. Tail 
particles extend to a radius of nearly 5 times the radius 
of the initial disc, despite the fact that the disturbance is 
impulsive. This value is in good accord with that predicted 
by equation (18) (as it should since both are realizations of 
the same approximation). The final panel shows the orbital 
excursions of groups of stars into the halo at a late time when 
the orbital phases at different initial radii have drifted too 
far to allow the formation of coherent structures. Although 
tail stars are perturbed to very radial orbits, they cannot 
form a recognizable tail on their second passage out into the 
halo. The lifetime of the tail is roughly the orbital period of 
the particles flung to the greatest distances. 

Similar results are seen in the milder case, A' = 0.5, 
shown in Figure 5. Here again the maximum tail length is 
well predicted by equation (18). The formation of azimuthal 
caustics on the leading edges of the tails is also visible in this 
figure. They form at about the time shown in the second 
panel, near the expected azimuths and 57r/4, see Sec. 

2.7). 

Caustic edges are even more visible in the still milder 
case, A' — 0.3, shown in Figure 6. Here the tails look more 
like spiral waves in an oval disc than tidal tails. All three 
of these figures have some general characteristics in com- 
mon. All are bilaterally symmetric as they must be given 
the symmetry of the simple tidal perturbation of equations 
(5) -(7). 

The tidal tails, and especially their compressed leading 
edges, are relatively straight in this approximation. Figures 

7 and 8 show that they are not as straight in the Keplerian 
potential. These figures also confirm the point emphasized 
by |Duc, Bournaud &c Masset| ( |2004| ), that with other pa- 
rameters equal, Keplerian tails are much longer than those 
of logarithmic potentials. Quantitatively, the maximum ra- 
dial extents are in good agreement with equation (18), with 

8 = 1/2. At some times the Keplerian tails look unrealis- 
tically broad compared to typical optical images. However, 
galaxies do not generally have rotation curves that decline 
as rapidly as the Keplerian. 

This effect is further illustrated in Figure 9, which shows 
the tails that develop after 5.0 time units in six model discs 
differing only in the exponent of their power-law potentials. 
(Compare the second to last panels in Figs. 6 and 7.) The 
upper left panel shows the case of the logarithmic potential 
(6 = 0), and the lower right panel show the Keplerian po- 
tential (5 = 0.5). The other panels have intermediate values 
of the exponent S and show a continuous increase in the tail 
length toward the Keplerian case. Otherwise the tails are 
qualitatively similar across the range of potentials. 
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Figure 7. Like Figure 4, but with a Keplerian (8 = 0.5) potential, 
and a tidal impulse amplitude of A' = 0.35 (at q = 1.0). Note 
that the appearance of the thin filaments in the middle panels is 
an artifact of the initial conditions. 




Figure 8. Same as Figure 7 (Keplerian potential), but with a 
tidal impulse amplitude of A' = 0.20 (at q = 1.0). 



Figure 9. Analytic model discs formed from a tidal perturbation 
of amplitude A' = 0.35 (at q = 1.0). Each panel shows a model 
with a different value of the potential exponent 8, which is given 
in the upper left corner. In all cases the model is shown at a time 
of 5.0 units. 



Figures 4-8 provide a large-scale view of the discs and 
tails, while we give one example of detailed structure within 
the disc in Figure 10. This figure shows results of a A' = 0.6 
amplitude encounter, with a potential intermediate between 
the logarithmic and the Keplerian (5 = 0.2). The model 
produces substantial tails. The inner disc is essentially an 
ocular structure (centred at about the azimuths predicted 
in Section 2.9), with strong internal spirals. Associated with 
the ocular and spirals the disc also produces strong swallow- 
tail (or higher order) caustics. E.g., paired ones with their 
strongest cusps located at about the 1 and 7 o'clock posi- 
tions in the figure (or about PA —30° and 150°). Presumably 
dissipative gas piles up in the converging flows that define 
the densest regions of the swallowtails, so these are good 



candidate sites for induced star formation (see e.g., Renaud, 
et al.|2009[|Renaud, Gieles fc Boily|20TT| . The outer caustic 
regions could be the locus of hinge clumps discussed below. 

Some caveats should be noted, however. First, the ex- 
treme symmetry of these structures is the effect of the as- 
sumed pure tidal perturbation, and is not realistic. Secondly, 
such structures are not common at this level of approxima- 
tion. Their production requires high amplitudes, very high 
amplitudes in the case of the logarithmic potential. These 
caveats are relaxed in the case of nonlinear, asymmetric per- 
turbations considered in Sec. 5. However, first we must con- 
sider another extension to the case of simple tidal perturba- 
tions. 
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Figure 10. Snapshot enlargement of a model disc with a poten- 
tial intermediate between the logarithmic and Kcplcrian (<5 = 0.2) 
at time t = 3.0 units, and an amplitude of A' = 0.60 (at q = 1.0). 
Paired, higher order caustics (e.g., swallowtails) arc evident in the 
outer disc near the base of the tails. 



3.2 Brief comparisons to numerical discs and tails 

Figures 1-3 and the discussion concerning them provide a 
basic view of the accuracy of the analytic models, but much 
more analysis would be required for a complete evaluation. 
That is beyond the scope of this paper, but following the 
theme of this section, we will pursue the issue somewhat 
further via another graphical comparison. This is provided 
by Figure 11, which was created for comparison to Figure 
5. All model parameter values are the same in both figures, 
including the tidal impulse amplitude. The difference is that 
following the impulsive perturbation to the initially circular 
orbits, the particle trajectories were computed numerically 
for Figure 11. Of course the result is not the same as a 
fully self-consistent, N-body simulation, but is a step in that 
direction, with self-consistent orbits calculated in the fixed 
potential. 

The results in Figure 11 clearly differ in some ways from 
those of Figure 5. The tails are somewhat longer at com- 
parable times (note the slightly different scales of the two 
figures), and longer lived in the numerical model. In terms 
of these properties, analytic models with larger tidal ampli- 
tudes (e.g., A' = 0.6) provide a better fit to the numerical 
results. However, given that the analytic model is formally 
accurate only to order e, the overall agreement is very good. 

Another clear difference is between the thicknesses of 
the tails in the two figures. A more precise description would 
be that the dense, but thin tails, in Figure 5 have not sep- 
arated from the parent disc, which is more elongated than 
that of Figure 11. We can understand this with the help of 
Figure 1, where in the 4 trajectories at the bottom of the fig- 
ure the analytic curves usually lie significantly above the nu- 
merical curves, i.e., at significantly larger radii. Nonetheless, 
these differences are again well within the expected limits of 
the approximation. 



4 EFFECTS OF TWO IMPULSES 

Calculations of two impulses separated by a finite time in- 
terval can provide the simplest model for a prolonged distur- 
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Figure 11. Same as Figure 5, except particle trajectories are 
computed numerically. See discussion in text. 



bance, and thus, the simplest representation of the difference 
between prograde and retrograde encounters between disc 
galaxies. In this section we will explore analytic approxima- 
tions to this case. Most of the relevant equations, and some 
limits are described in the appendix. 



4.1 Prograde versus retrograde interactions 

The equations derived in the appendix show that, in op- 
timal circumstances, the apoapse radius of a stellar orbit 
can grow geometrically with successive impulses. What is 
required for this optimal amplification? Ultimately it is due 
to two distinct physical effects, which correspond to the nu- 
merator and denominator of the term in square brackets in 
equation (45). The first effect, codified in the numerator of 
the term in equation (45), is the increase in angular mo- 
mentum resulting from the prograde perturbation (see eqs. 
(6) and (38)). More specifically, it concerns the region where 
stars receive the maximum increase in angular momentum. 
In the first impulse this was along the azimuthal line called 
the backbone in Section 2.5. Clearly, if stars on the back- 
bone also receive the largest angular momentum increase in 
the second impulse their amplitude and apoapse radius will 
be relatively large. Such stars must lie on both the first or 
second backbones, i.e., at azimuths of — n/2 or 37r/4 relative 
to the closest approach point at both impacts. 

This is most easily achieved if these stars rotate at the 
same average angular frequency as the companion (or equiv- 
alent angular frequency in the case of hits by two different 
companions). Since angular frequency in the disc generally 
varies with radius, generally at most one spot of the origi- 
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nal backbone (call it a particular vertebra) will satisfy this 
condition, and be found on both backbones. 

For most stars, their relative azimuth at the second im- 
pact will be different from the first. They will not get the 
maximum boost from both, and in some cases the angular 
momentum impulses may be of opposite signs. This is the 
interference between the two impulses noted above. In the 
case of a retrograde interaction, it is practically impossible 
for a star on the original backbone to find itself on the sec- 
ond. (Unless, of course, the companion and the disc stars 
have combined rotations of more than 2n radians between 
the two impulses.) The retrograde case also generates an in- 
terference pattern of velocity impulses on the disc, but one in 
which the original winners (in terms of angular momentum 
gains), are much less likely to win again. 

The second factor, for determining the amplitude of a 
stellar orbit after the second impulse, is given by the denom- 
inator of the bracketed term of equation (45). Physically, the 
value of this term is determined by the phase along the p- 
ellipse orbit at the time of the second impact. I.e., if the star 
is at apoapse it will be most strongly affected. The optimal 
effect, described by equations (47) and (51), occurs for stars 
that lie on both backbones and happen to be at apoapse at 
the time of the second impact. 

To get the maximal input from this second factor, the 
azimuthal phases of the p-ellipse immediately after the first 
impulse and just before the second, must be given by, 



6i2 = On + cA(f) : 



(29) 



Using equation (12) we solve for the total azimuthal phase 
change, 



A(f> : 



TV 

m 



(2(j> + | + V) • (Modulo 2tt.) (30) 

The two amplitude effects are not directly connected, so 
such a conjunction between them would appear to be a co- 
incidence. However, there is a bit more to it than that. The 
second effect requires a significant time between impulses, 
while the first effect requires co-rotation between a region 
on the backbone and the companion. 

If the two impulses represent a prolonged encounter, 
then the coincidence of the two effects does not seem un- 
likely. At the least, the continued angular momentum in- 
put into a particular region on the initial backbone, lasting 
through half an orbital period of a kind of average p-ellipse 
orbit in that region, seems reasonably likely. It is certainly in 
accord with many published simulations which show that the 
longest tails are produced by prolonged prograde encounters, 
and is the simplest analytic representations of such encoun- 
ters (see also D'Onghia, et al.|2010 and Nelson & Tremaine 
1999). The present treatment shows clearly the effects of 



resonance in the azimuthal frequency. 



4.2 Conjectures on mechanisms of tidal dwarf 
formation 

It is likely that there are several specific mechanisms for 
producing dwarf galaxies in tidal tails. One of these may 
be based on the azimuthal pileup discussed in Section 2. 
However, a number of the more substantial tidal dwarfs seem 



to lie at or near the end of tails. This requires a mechanism of 
radial pileup. Nothing like that was evident in the impulsive 
tails described in the previous sections. On the other hand, 
the phenomenon described in the previous subsection of a 
'resonant vertebra' along the initial backbone of the tail, 
which receives the strongest azimuthal velocity input in the 
disc from a prolonged encounter, could lead to such an effect. 

With greater angular momentum than surrounding re- 
gions, stars in the resonant vertebra region could overtake 
others in the tail, which were initially located at larger radii, 
leading to a density enhancement. In the case of gas clouds 
this could lead to an actual pileup or shock wave. Both ef- 
fects could lead to local gravitational instability and the 
formation of a self-gravitating tidal dwarf. Regions interior 
to the resonant vertebra would be tidally stretched more 
strongly than in single impulse tails, and as a result, would 
be less likely to produce large star clusters. However, there 
could be a competition with the effects of azimuthal pileup, 
though this is likely to occur at an earlier time than the ra- 
dial pileup. These considerations could explain the observed 
tendency for tails with relatively large dwarfs at their ends to 



have few clusters along their length, and vice versa (Knier- 



man, et al.|2003| |Mullan, et al.||2011| |Kaviraj, et al.|20~12 i 

They could also help explain the rarity of large dwarfs 
at the end of tails. If the resonant vertebra is located at the 
outer edge of the initial disc (or beyond), there will be no 
strong pileup. If the resonant vertebra is located deep within 
the initial disc, then the pileup will more likely occur within 
the disc rather than well out in the tail. This might instead 
play a role in the phenomenon of hinge clumps (as defined 
in|Hancock, et al.|[2009| |Smith, et al~1|2010a[ |Smith, et al~1 
|2010b[ >. 

In any case, the discussion in these last two subsec- 
tions illustrates the exploratory power of the two-impulse, 
p-ellipse model for disturbed stellar orbits, and their collec- 
tive behavior. 



4.3 Graphical realizations with multiple hits 

Two sample realizations of tail evolution with multiple hits 
are shown in Figures 12 and 13. In both these models there 
are three successive (prograde) tidal impulses, each with an 
amplitude of A — 0.35. The formalism for two hits given in 
the preceding sections is straightforwardly extended for any 
number of impulses, and readily implemented in a computer 
script. The choice of three impulses here is arbitrary. It is 
sufficient to illustrate the effect of multiple hits and compare 
to the formulae above, but not as complex as a simulation 
of a very prolonged encounter with numerous impulses. 

In the case shown in Figure 12 the effective azimuthal 
frequency of the companion relative to that of a disc star is 
about 3.7q. This ratio is equal to one at an initial disc radius 
of q = 0.27. The resonant vertebrae of the backbone should 
be located at about that radius, rather deep inside the disc. 
The azimuths of the three impacts are 0°, 70° and 140°. 

For the case shown in Figure 13 the effective azimuthal 
frequency of the companion relative to a disc star is q. 
Clearly, this ratio is equal to unity at a disc radius of q = 1.0, 
at the outer edge of the disc. In this case, the azimuths of 
the three impacts are 0°,57° and 114°. 

Equation (47) suggests that the maximum amplitude in 
these cases is about 0.9. The approximate equation (51) sug- 
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Figure 12. Snapshots of the evolution of a disc experiencing 
three tidal impulses, each of amplitude A' = 0.35 (at q = 1.0). 
The impulses occur at azimuths of about 0°,70° and 140°, and 
at times of 0.05, 0.75 and 1.5 units. 



gests that the maximum obtainable apoapse radius is about 
1.7 2 = 2.9 times the initial radius. The better approximation 
of equation (50) suggests that this factor is more like 4.0. In 
Figure 12 the tail does not stretch out to such a large factor 
relative to the initial disc radius. The maximum radius in 
this case is somewhat larger than 3. This is not surprising 
given that the strongest orbital resonances are with stars lo- 
cated well within the initial disc. In the third panel of Figure 
13 the tail does stretch out to a radius of 4.3, which might 
be expected since the orbital resonance is with stars at the 
outer edge of the disc. 

At their longest extents the tails in both Figures 12 and 
13 are more like those in Figure 4 (single impulse amplitude 
of A = 0.8) than those of Figure 6 (single impulse ampli- 
tude of A = 0.3), which emphasizes the cumulative effect of 
the multiple impacts. However, like low amplitude single im- 
pulse cases, the tails of Figure 13 are smooth, with no major 
clumps. The fourth panel of Figure 12 shows clumps in the 
form of swallowtail caustics in both tails (e.g., at time 3). 
This appears to be an example of the mechanism for tidal 
dwarf formation discussed in the previous subsection. 

Given the number of free parameters in these models 
(e.g., impact amplitudes, spacing and frequency), or those 
in more realistic prolonged encounters, it is clear that a very 
wide variety of tail structures could be produced. Further 
exploration of these forms is beyond the scope of this paper. 
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Figure 13. Same as Figure 10, except that the impulses occur 
at azimuths of about 0°,57° and 114°, and at times of 0.05, 1.0 
and 2.0 units. 



5 NONLINEAR TIDAL DISTURBANCES 

Rather than trying to explore the evolution of ensembles 
of multiply perturbed stellar orbits, we want to return to 
some issues remaining from the discussion in Sections 2 and 
3, specifically, the unrealistic symmetry of the tidal tails il- 
lustrated there. Because p-ellipses have been shown to give 
good approximations to orbits in the power-law potentials 
considered here, we do not believe that the fault lies with 
these orbit functions. It appears that we need a more accu- 
rate representation of the collision impulses than provided 
by equations (5) and (6). These equations are only linear 
in the factor q/r m i n , and so are not likely to be sufficiently 
accurate for either strong or prolonged encounters, which 
will usually be close encounters in both cases. More com- 
plex versions of the tidal impulses will give more complex 
expressions on the right hand side of equation (10). However, 
this does not pose a great difficulty. 



5.1 Nonlinear Analytic Models 

To consider the effects of nonlinear perturbations we will use 



the functions of Gerber & Lamb ( 1994 ) for the impulsive 



velocity disturbance from a Plummer potential in place of 
equations (5) and (6) above. ( |D'Onghia, et al.| ( |2010[ > give 
more accurate representations for various cases. However, 
the former will suffice for the present purpose of illustrating 
the effect of nonlinear terms in the ratio of the initial stellar 
radius to the closest approach distance.) Specifically, we use 
the velocity impulses given by equations (7a) and (7b) of 
Gerber and Lamb, with softening length 7 = 0, which can 
be written as, 
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Av r = Aryv^o 



cos (2(j> ) — 77 cos (4> ) 
1 — 2r/cos (<f> ) + rj 2 



sin (20 o ) — 77 sin (<j> ) 
1 — 2?7cos (<j> ) + rf 2 



Aw^ = ArjVcfto 
where for convenience we adopt the variable 



(31) 



(32) 



(33) 



With the substitution of these expressions, equation (10) 
becomes, 



( \ + 8) me sin (m0 o ) 



[1 + ecos (m0 o )] 

A' (cos (2<p a ) — rj cos (<j> )) 
1 — 2??cos (<j>o) +f]' 2 — A' (sin (2<j> ) — 77 sin (<p )) 



(34) 



and the orbital solutions are still of the same form as equa- 
tions (11) (13). That is, equation (11) is unchanged, but 
equations (12) and (13) for e and ip become, 




5 10 
Azimuthal Advance 



Figure 14. Same as Fig. 1, but with a nonlinear, impulsive dis- 
turbance as given by equations (31) (33) with rj = 0.10 when 
g = 1.0. 



^ + S^j me = [ A' 2 (cos (2<p ) - 7700s (<p )) 2 + 

(2?7cos(0 o ) - V 2 + A'(sin(20 o ) - rj sin (<f> ))) 2 ]s , (35) 
and, 

tan (20 o + %l>) = [ \ + S) mx 



(2r?< 



b ) - -rj 2 + A' (sin (2<f> ) - r;sin (0 O ))) 



A' (cos(20 o ) - 7j COS (<j>o)) 



(36) 



The expression for the eccentricity in equation (35) is always 
positive, but can exceed unity. 

Equation (35) is a generalization of equation (12). Re- 
call that orbits with eccentricities given by equation (12) 
were found to be moderately accurate (see Figures 1 and 
2). Subsequently, equation (12) was modified with the am- 
plitude dependent coefficient to the cosine term shown in 
equation (15) to moderate unrealistically large eccentricities 
in large amplitude orbits. A second modification applied in 
going from equation (15) to (16) was used to increase the 
unrealistically low eccentricities of orbits near the backbone 
(Figure 3) . Given the complexity of the perturbation (equa- 
tions (31) and (32)) it is not surprising that the orbital solu- 
tions derived from equation (35) are somewhat less accurate 
over the whole range of A' and 8 values considered than in 
the pure tidal case. Adding the extra coefficients of equation 
(16) to equation (35), is also not as effective in correcting 
the orbits as it was in the tidal case. However, we can do 
nearly as well by using slightly modified coefficient. By ex- 
perimentation we have found the following expression for the 
eccentricity works well over the relevant parameter range, 



+ 8 me - 



- + 5} m 2 x 



A" I I y^Ta? ) cos (20 o ) - r? cos (0 O ) 1 +U + 



[2rj cos (<p ) — t 2 + A' (63 sin {2<j> ) — rj sin (<p ))] 2 ] 



with 63 = (l + A') 



1+6S- 1 



(37) 



Figure 14 shows a comparison between the resulting analytic 
orbits and numerically integrated orbits. 



5.2 Results from the nonlinear model 

The first result from the nonlinear model of the last subsec- 
tion is a symmetry breaking. In the basic model of Section 2 
the amplitude of the tidal perturbation was determined only 
by the factor A' , in the nonlinear model there is an addi- 
tional dependence on the factor rj (see equations (31) - (33)). 
The nonlinear effects are most important for stars with large 
values of rj, i.e., for small impact parameters and for stars 
with relatively large initial radii q. In the limit of small rj, 
the nonlinear effects are negligible, and the orbit function 
is nearly the pure tidal p-ellipse. In this subsection we will 
look at several example stellar discs, like those discussed in 
Section 3 for the tidal case. 

Figure 15 shows a series of snapshots for a case similar 
to that shown in Figure 5, and with rj — 0.1 for stars on the 
outer edge of the disc. For example, this could represent a 
case of a companion galaxy flying rapidly past (impulsive 
collision) a 6 times more massive primary galaxy with a 
closest approach distance of 10 times the radius of its stellar 
disc. Despite this relatively large distance of closest approach 
(in terms of the radius of the small galaxy), the asymmetries 
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Figure 15. Snapshots of a model in the logarithmic potential 
with an amplitude of A' = 0.60 (at q = 1.0), but with a nonlinear, 
impulsive disturbance as given by equations (31) (33) with 7] = 
0.10 when q = l.Q (i.e., closest approach at 10 times the initial 
disc radius). 



generated by the nonlinear perturbation are obvious. The 
general size and structure of the two tails are not greatly 
different. However, the caustic waves are more complex in 
the nonlinear cases than in the tidal cases. By the last times 
shown the asymmetries have developed into a one-armed 
form. 

Figure 16 shows a lower amplitude case, with A' = 0.3 
(at q — 1). The closest approach distance is also smaller 
(rj = 1/3), so the nonlinear asymmetries are larger in this 
case. The bilateral symmetry of the pure tidal part of the dis- 
turbance is hardly evident. (Compare to Figure 6.) A series 
of different caustic waveforms appear (as in some multi-hit 
cases of Section 4) within the one-armed wave. Interestingly, 
at some times (e.g., times 3 and 5 in the figure) there are 
over-dense regions at the end of the arm. Evidently such a 
seed for the growth of a tidal dwarf galaxy can be quite read- 
ily produced in an impulsive, but nonlinear, encounter. The 
over-density is produced by caustic waves, but via a mech- 
anism that is different than that conjectured in Section 4.2. 
Note that the mechanism requires a two-dimensional inter- 
action, otherwise the crossing star-streams could be in dif- 
ferent planes, and not contribute to a true density enhance- 
ment. 

High-density zones of converging flows are also common 
at the base of the tails in Figures 15 and 16. These provide 
models of a hinge clump. Like many observed hinge clumps 
(see Smith, et al.|2010b l, and other features in the nonlinear 
models these are not generally symmetric. In tidal models 
they do have a twin on the opposite side of the disc (as 
seen in Figure 4 for example). Hinge clumps are a recently 



Figure 16. Same as Figure 13, but with A' = 0.6 and r\ = 0.333 
(at q = 1.0). 



identified structure, with as yet only a few known examples. 
The existence of a production mechanism, like that shown 
in the figure, provides support for the idea that they are a 
true physical structure, and not simply random clumps of 
star formation. We will consider more comparisons to obser- 
vation in the next section. 

Figure 17 shows a moderate disturbance (A' — 0.3), 
with a moderate distance of closest approach (rj = 1/5), in 
the Keplerian potential. Asymmetries like those in Figures 
15 and 16 are apparent. The most impressive feature is the 
long, one-sided tail that develops at late times. 



6 COMPARISONS TO OBSERVED FORMS 

In developing the models presented in this paper we have two 
general goals: 1) to provide simple analytic tools for study- 
ing the theory of tidal tails and associated phenomenon in 
general, and 2) to provide analytic tools to aid in deciphering 
the collision history of individual observed systems. The pre- 
vious sections have been occupied with the first goal. In this 
section we focus on the latter, and specifically, undertake the 
relatively modest goal of looking for systems that the various 
models might be applied to, and especially some tentative 
examples that the single-hit, impulsive models might apply 
to. 

To begin we specify candidate impulsive cases as those 
with a visible companion that is well separated, and exclude 
systems that appear to be in prolonged interactions or im- 
minent mergers. Among the former, symmetric cases may 
be more purely tidal, asymmetric ones nonlinear. For candi- 
dates chosen within these restrictions, the following relations 
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Figure 17. Same as Figure 14, but in the Keplerian (8 = 0.5) 
potential, and with r] = 0.2 (at q = 1.0). 



predicted by the models can be examined. 1) The tidal dom- 
inant cases should have tails that are generally smooth along 
their length. 2) They may show evidence of a leading edge 
caustic. 3) In all cases long tails should be associated with 
large companions, especially in the tidal case where they of- 
ten do not approach as closely as in the nonlinear case. 4) 
The longest tails have the least curvature. (This is true even 
in the Kepler potential.) 5) Irregularities within tails, pro- 
duced in nonlinear interactions, should be associated with 
other asymmetries in the disc. 

We note that observed tail characteristics depend on 
the observational waveband, and whether the observed con- 
stituent consists of gas, young stars, or old stars. However, 
the properties described above are very basic geometric ob- 
servables that relate quite directly to model characteristics. 
We also note that the available and adequately resolved ob- 
servations are not sufficiently extensive to make detailed 
statistical comparisons (though global statistical properties 
have been studied recently by|Smith, et al.||2007| |Smith, et| 



al.|2010a||Mohamed, Reshetnikov fc Sotnikova||2011[ ). 

With these caveats in mind we will consider example 



systems that are chosen from the Arp Atlas (Arp 1966), 



which provides an optically well-resolved sample of sufficient 
size to contain examples of systems that are somewhat rare 
even among the rare class of strongly interacting galaxies. 
The examples discussed below were selected by eye, initially 
independently by the two authors, with subsequent consul- 
tation on systems not chosen by both. This is not a rigorous, 
objective procedure, but it accomplishes the goal of finding 
tentative examples and allowing simple tests of the relation- 
ships suggested by the models. 



6.1 Long, straight tails 

One interesting aspect of the impulsive, pure tidal models 
shown in Figures 4, 5, and to some degree in Figure 7, is 
the relative straightness of the long tails produced in those 
cases. (At least in the parts of the tails that lie some dis- 
tance from the disc.) We expect them to be rare, because 
a strong perturbation is required to produce a long tail, 
but an impulsive encounter is needed to produce the sim- 
ple, symmetric form. In fact, long, symmetric tails are very 
unlikely, because a strong perturbation generally requires a 
large companion, and at most a slightly hyperbolic relative 
velocity. In such cases, dynamical friction can hold the col- 
lision partners together, with the formation of a bridge plus 
counter-tail, rather than two symmetric tails. 

These expectations are borne out by the recent statisti- 
cal study of the lengths and curvature of tidal tails by |Mo-| 
hamed, Reshetnikov & Sotnikova (20111. They found that 



most tidal tails have lengths less than three times the di- 
ameter of the disc; only a small fraction have lengths longer 
than four disc diameters. Furthermore, most tails have sig- 
nificant curvature on the sky, thus intrinsically straight tails 
are likely relatively rare. 

When observed (doubly or singly), long, straight tails 
are often assumed to be edge-on. For example, the straight- 
ness of the northern tail of Arp 242 (the Mice) may be caused 
by the disc being edge-on (e.g., Toomre fc Toomre]|1972| 



Mihos, Bothun fc Richstone 1993 Bournaud, et al.||2004 



In general, however, the tendency of tails to stay in their 
initial plane means this is unlikely unless the parent disc is 
also observed edge-on. (Not to mention the fact that it is 
statistically unlikely.) 



6.1.1 Symmetric objects 

Nonetheless, the Arp Atlas does have a few long tailed sys- 
tems that are close to symmetric, including objects: 64, 65, 
99, 295, and 296 (see the top row of Figure 18). The disc 
of Arp 295 does look edge-on, while that of Arp 296 is too 
small to tell. The discs of the first three objects do not ap- 
pear edge-on. In all of these cases, except Arp 64, the tailed 
galaxy has a massive companion, at a distance of at least a 
couple of disc diameters, consistent with expectations from 
an impulsive, tidal model. Arp 64 has a hook at the end 
of one tail, and may have two companions, and so, may be 
more complex dynamically. In Arp 64, Arp 65 and Arp 99, 
the tails connect to strong symmetric spirals within the disc. 
The small fraction of these few long symmetric tails out of 
the few hundred systems in the Arp Atlas, seems qualita- 
tively consistent with the restricted class of encounters that 
can produce them. 



6.1.2 Asymmetric objects 

We can also find apparent examples of long, straight, but 
single or very asymmetric tails, including Arp objects: 100, 
101, 102, and possibly 190, if it is not presently merging (see 
the bottom row of Figure 18). These tails suggest a very 
nonlinear disturbance with at least one strong impulse like 
the models shown in Figure 16 and 17, but on a larger scale, 
and at somewhat later time. It is true that such singular tails 
are easier to manufacture in the Kepler potential than in the 
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Figure 18. Top row: Examples of Arp galaxies with long straight symmetric tails. Left to right: Arp 64, 65, and 99. Bottom row: 
Examples of Arp galaxies with long, straight, but asymmetric tails. Left to right: Arp 100, 101, and 102. Images from |Arp| ( |1966p 



logarithmic, so in the Arp galaxies they may extend beyond 
the point where the parent galaxy halo potential begins to 
fall off. 

It is unfortunate that these particular examples are ap- 
parently quite distant, so it is hard to discern caustic struc- 
tures within their discs. The tails of both Arp 100 and 101 
appear smooth, without prominent clumps, yet they are nar- 
row, with blue optical colors and high optical surface bright- 
nesses ( |Schombert, Wallin fc Struck-Marcell|1990 |. 



In contrast, the ends of the long straight tails of Arp 
102 and 105 are surrounded by a lot of debris, with knots 



of star formation at the ends (Schombert, Wallin & Struck- 



Marcell|1990l|Duc, et al.|1997| >. The base of the Arp 105 tail 
is narrow and well-defined, while Arp 102 has a broad tail. 



Due, et al. ( 1997 1 suggest that the straightness of the Arp 



105 northern tail was caused by the interaction being seen 
almost edge-on. The Arp 190 tail has star formation near 



its tip (Schombert, Wallin & Struck-Marcell 19901 and a 



broader base. Arp 179 does not have an obvious companion 
galaxy, thus it may be a merger remnant. Many of the long 
straight tails in the Arp Atlas, like those of the Arp 226, the 
Atoms for Peace galaxy, are in systems where there is clear 
evidence that a merger is in progress (e.g., |Toomre|p78| . 

A final cautionary example is provided by Arp 188 (the 
Tadpole), which has a long straight tail. However, Hubble 
Space Telescope (ACS) imagery seems to affirm the impres- 
sion from earlier images that the disc of this galaxy has 



been distorted in three dimensions (Tran, et al. 2003 1. Extra- 



planar disturbances and their consequences have not been 
considered in this paper. In cases when such disturbances 



are small the model formalism could easily be extended to 
include epicyclic motions in the vertical direction. 

The fact that we have to look out over large volumes to 
find examples of this class is suggestive. It must be harder 
to produce nonlinear disturbances impulsively, than pure 
tidal disturbances. Nonlinear disturbances will usually be 
the result of prolonged encounters. 



6.2 Long, somewhat curved tails 

If we relax the straightness requirement, and allow some- 
what more curved tails, we find more systems, including 
Arp 33, Arp 173 and Arp 174 (see top row of Figure 19). In 
the case of the first of these most of the curvature is near 
the base of the tail. In going from Figure 4 to 6 we see a 
trend of more curvature at the base of the tail with decreas- 
ing perturbation amplitude. The companion in the Arp 33 
system appears smaller than the long-tailed object, suggest- 
ing a smaller perturbation than in the systems considered 
above. In the SDSS images the long-tailed galaxy is quite 
blue, with prominent knots of star formation near the end. 
The companion also has a nice, symmetric pair of shorter 
tails, but without bright star formation regions. 

The tailed galaxies in the Arp 173 and 174 systems are 
both close to larger companions (in projection), so it seems 
likely that the encounter has been prolonged. In contrast to 
the Arp 33 tails, the Arp 173 tail is smooth, with relatively 



red UV/optical colors ( |Smith, et aL]|2010a[ ). The Arp 174 
tail is relatively diffuse, also without prominent star forming 
regions. 

Another system with long somewhat curved tails is Arp 
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Figure 19. Top row: Examples of Arp galaxies with long but somewhat curved tails. Left to right: Arp 33, 173, and 174. Bottom row: 
Examples of Arp galaxies with long very curved tails. Left to right: Arp 58, 72, and 82. |Arp| ( |1966^ 




Figure 20. Examples of Arp galaxies with short symmetric tails. Left to right: Arp 31, 285, 286. Arp (1966) 



96. The tails are faint and diffuse, while the inner spiral 
arms are quite prominent and curve back on themselves. The 
tails appear to branch out from the spiral arms at about the 
location where the arms start curling back inwards. A bright 
elliptical lies near the end of one of the tails. 



6.3 Long curved tails 

The Arp atlas contains a number of long very curved tails, 
many of whose parent galaxies seem to be near merging 
or merged. There are also a number whose companions are 
more than a couple diameters distant, so merger does not 
seem imminent. These include Arp objects: 58, 72, 82, 97, 98, 
252, and possibly Arp 35, though its disc looks sufficiently 
disturbed that it may have experienced a recent merger (see 



bottom row of Figure 19). The results of many numerical 



merger models at pre-merger times (e.g., Barnes & Hern- 



quist|1992b1|Mihos, Dubinski fc Hernquist|1998[|Duc, Bo~ 
naud & Masset 2004|, and a numerical model for the Arp 
82 system in particular (see Hancock, et al.|2007 l, teach us 
that such tails can be formed in prolonged interactions with 
a bound companion. It takes a considerable time to pull tails 
out to the observed lengths, yet the companions in these sys- 
tems are relatively nearby. This is circumstantial evidence 
for the conclusion that the companions are bound, and have 
been interacting for some time. These types of tails are not 
produced in the single impulse models described above. In 
the multiple hit cases shown in Figures 12 and 13 the inter- 
action is probably not strong enough, or prolonged enough, 
to produce them. 



Tidal tails in galaxies 19 

















[a 











Figure 21. Some example of hinge clumps, as seen in the optical 
and the UV. Top row: The western galaxy of Arp 65. Bottom row: 
The southern galaxy of Arp 305. Left column: Optical images 
from the Arp Atlas. Right column: GALEX NUV images from 
|Smith, et aT] ( |2010a| |. North is up and east to the left. In both 
systems, two tails extend to the northwest and southeast. In the 
Arp Atlas pictures, note the 'loop-back' segments from the tails 
curving back into the inner regions. Also note the two UV-bright 
'hinge clumps' where the 'loop-backs' diverge from the tail. 



6.4 Shorter symmetric or nearly symmetric 
systems 

Long tails are the most spectacular tidal remnants. Short 
tails are also of interest, though it is not always entirely clear 
that they are tidal products. It is possible that some may 
be spiral waves in relatively low surface brightness discs. On 
the other hand, strong spiral waves seem unlikely in a low 
surface density disc, unless there is some external driver. In 
systems for which there is little data available beyond atlas 
images ambiguities remain. 

The Arp atlas provides several examples of short, but 
quite symmetric tails, i.e., in objects: 1, 9, 31, 77, 80, 83, 
84, 214, 285, 286, 293 and 316 (see Figure 20). All of these 
galaxies have possible companions, most with roughly the 
same redshift, though the apparent mass ratios extend over 
a large range. Four of these objects (Arp 9, Arp 77, Arp 
80 and Arp 285) have obvious bars, so that the 'tails' could 
possibly be bar driven spirals. Five of these objects, and 
among the most symmetrical, are themselves the smaller 
companion. These objects are: Arp 83, Arp 84, Arp 286, 
Arp 293, and Arp 316. 

We have previously published numerical models of two 
of the systems above, Arp 84 (Kaufman, et al. 1999 also 



see Puerari, Valdez-Gutierrez & Hernandez-Lopez 2005 \ and 
Arp 285 ( |Smith, et al.|2008[ ). In both cases the models sug- 
gest that the galaxies have experienced a fairly rapid (im- 
pulsive) flyby, and are now moving out to apoapse. In the 
case of Arp 84 the models of Kaufman et al. suggest a closest 
approach near the outer edge of the primary disc, on a tra- 



jectory that is slightly tilted with respect to the plane of the 
primary disc. The inclination of the companion disc relative 
to the primary galaxy is similarly modest. The encounter in 
the simulation was quite impulsive. In the simulation for Arp 
285 the relative orbit is in the plane of the primary (tailed) 
galaxy disc, and the flyby is rapid, so the simulation is very 
much like the models above. 

Many of these systems bear strong resemblances to the 
middle panels of Figures 5, and 6, and the timescales of those 
panels are also about right (e.g., based on the numerical 
models). Stated differently, the analytic models imply that 
short, symmetric tails appear shortly after impulsive, but 
dominantly tidal encounters, and these observational and 
numerical examples are in accord with those correlations. 



6.5 Swallowtails and hinge clumps 

As discussed in Section 3 and as illustrated in Figure 10, 
the models sometimes produce swallowtail caustics near the 
base of tidal tails. Intersecting streams of gas may trigger 
star formation at these locations, perhaps accounting for 
the 'hinge clumps' observed in some interacting systems. 
Two galaxies that match the basic morphology of Figure 10 
are shown in Figure 21: the western galaxy in Arp 65 (top 
row) and the southern galaxy in Arp 305 (bottom). The left 
panels show expanded views of the Arp Atlas images of these 
two galaxies, while the panels on the right display matching 
GALEX NUV images. In the Arp Atlas image, as discussed 
by Hancock, et al. (20091, the Arp 305 galaxy has an ocular 



structure, with two fairly straight tails coming out of the 
ends of the 'eye'. Partway along the tails, luminous knots 
of star formation are seen ('hinge clumps') on the GALEX 
image. Near the hinge clump in the southern tail, an optical 
loop is seen in the Arp Atlas picture that arcs back into 
the main galaxy. A similar, but less prominent, loop-back is 
also seen near the northern hinge clump. This morphology 
resembles the model in Figure 10. 

A similar structure is seen in the western galaxy of Arp 
65, except that the inner ocular is not completely closed, 
being an open symmetric spiral. The two hinge clumps, seen 
to the northeast and southwest of the nucleus, lie near the 
location where 'loop-backs' separate from the tail in the Arp 
Atlas picture. These hinge clumps are also seen in the Spitzer 
8 micron images of |Smith, etH] ( [20071 ). 

Another system that may have a similar structure is Arp 
99, which also has long straight tails like Arp 65. A published 
Ha map shows numerous knots of star formation along the 
tails, with particularly luminous clumps at the bases of the 
two tails (Igle sias-Paramo fc Vilchez||200l| ). The Arp Atlas 
image of this galaxy also suggests 'loop-backs' near these 
hinge clumps. Hinge clumps are also seen in galaxies with 
long curved tails. For example, Arp 82 also has a luminous 



hinge clump at the base of its northern tail ( Hancock, et al 



2009 1. Arp 58 may have a similar structure, judging from the 
Arp Atlas photograph and an archival GALEX UV image. 
Although more disturbed, UV and optical images of Arp 72 



(Smith, et al. 2010a Arp 19661, and infrared and optical 
images of Arp 240 ( |Smith, et al.||2007| l, suggest the same 
basic structure. 
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6.6 Observational comparison summary 

About 5-10% of the objects in the Arp atlas may provide 
examples of interacting systems that can be described by 
the impulsive analytic models above. This percentage might 
be considerably increased when the models for multiple or 
prolonged flyby encounters are included. This is about what 
we would expect given that the models of this paper are 
restricted to planar, flyby collisions, and exclude mergers (in 
progress or completed), and three-dimensional encounters. 

The sample systems also seem to be in general agree- 
ment with correlations suggested by the analytic models, 
noted at the beginning of this section, and discussed in the 
preceding subsections. We conclude by mentioning a cou- 
ple more. Firstly, we emphasize that most of the long tails 
of Section 6.1 and the short tails of Section 6.4 are quite 
smooth along their length, with the exception of some knots 
of star formation. The barred objects discussed in the pre- 
vious subsection are somewhat less smooth than the un- 
barred systems. This makes sense in the context of inter- 
ference effects between tidal and bar-driven perturbations. 
The shorter tails tend to have higher surface brightnesses, 
and some appear to have quite sharp leading edges, as we 
would expect for azimuthal caustics (e.g., Arp 83, Arp 285, 
Arp 293, and perhaps Arp 316). We would reemphasize that 
among the symmetric tails, the shorter ones have more cur- 
vature, and the long ones tend to be straighter. In sum, the 
models provide fair representations of the basic phenomenol- 
ogy, and the observations of those systems do not obviously 
contradict model properties. However, the comparisons here 
are only qualitative, and given the limited data it is hard to 
be sure that all the examples satisfy the approximations on 
which the models are based. 



7 SUMMARY AND CONCLUSIONS 

The p-ellipse functions describe orbits of stars in power law 
potentials very naturally. Because they allow large radial 
excursions, tidal tails and related structures are also natu- 
rally produced with these orbit functions. These model tidal 
structures appear to match the observations well in systems 
where the model approximations (e.g., of an impulsive pla- 
nar flyby) may obtain. We conclude with a summary of our 
main results. 

1. The p-ellipses provide simple, yet accurate approxi- 
mations to impulsively perturbed orbits. However, at mod- 
erate to large eccentricities the best fit is achieved with a 
change in the p-ellipse eccentricity from that in the simplest 
solution of equation (12), to the modified version of equation 
(16). 

2. One aspect of this fit, the result from equation (12) 
that e fa A', is an especially simple result, if not directly 
observable. 

3. The form of the tidal perturbation suggests that tails 
form around a backbone at a specific azimuth. The p-ellipse 
orbits for stars on the backbone provide a very simple model 
for the evolution of the simplest tails. This analysis yields 
a formula for how the maximum length of the tail depends 
on the amplitude of the tidal perturbation. This result could 
be estimated from angular momentum conservation without 
reference to the p-ellipse. However, the dependence on the 



halo potential is a simple result of using the p-ellipse equa- 
tion. Moreover, as shown in Section 3 the p-ellipse orbits 
can also reveal the structure of the whole tail as a function 
of initial position and time. (Also see Lynden-Bcll 2010b 



The generic evolution of the orbit ensembles also 
supports some observable generalizations (and accords with 
some simulation results). For example, tails produced by 
tidal impulses are smooth, symmetric, but generally develop 
an azimuthal caustic region of enhanced density on their 
leading edge, after traveling 90° from their initial azimuth. 
Long tails are very straight. Nonlinear (not purely tidal) 
impulses can produce higher order caustics, and less regular 
tails, as can multiple tidal impulses. In both cases, detailed 
tail structure depends on the halo potential (assumed rigid 
in these analytic models). 

5. The p-ellipse approximation provides the basis for 
interpretations of a variety of other tidal phenomena, in- 
cluding ocular rings (radial caustics corresponding to the 
angular caustics in the tails), and hinge clumps (as higher 
order caustics located near the base of tails). The orbit for- 
mulae can also be used to derive evolutions and scalings for 
these structures via analysis of the caustics. 

6. In cases of two or more impulsive encounters, the sin- 
gle impulse results can be extended to produce estimates of 
tail length in prograde versus retrograde cases, as a func- 
tion of the form of the potential. Two impulsive encoun- 
ters (e.g., viewed as a simple approximation to an encounter 
of finite duration) provide the simplest model for the pro- 
grade/retrograde dichotomy. 

7. In the case of two encounters, the map of the sec- 
ond perturbation across the disc interacts with the evolved 
velocities from the first perturbation in a way that can be 
viewed as an interference pattern. On the (second) backbone 
an especially strong perturbation can result from this inter- 
ference at a specific location. This could provide the seed for 
the formation of a tidal dwarf galaxy, as discussed in Section 
4.2. Further work is needed on this conjecture. Nonlinear or 
prolonged interactions can also generate high order caustic 
density enhancements that could seed tidal dwarf formation. 

8. The most general case considered in this paper is that 
of multiple nonlinear impulses (see appendix). We have not 
investigated this case, but the discussion in the appendix 
provides a template for studying it with any nonlinear per- 
turbation function. 

More generally, the p-ellipse orbits provide a useful tool 
for furthering our understanding of a wide range of phenom- 
ena in colliding or otherwise disturbed galaxies. Previous 
studies have considered the properties of the p-ellipses, and 
their unusual accuracy in approximating galaxy orbits (S06, 
Lynden-Bell| |2010a[ ) . The applications above help demon- 
strate how the family of functions may prove useful in many 
problems in galaxy dynamics. 
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8 APPENDIX: EQUATIONS FOR THE CASE 
OF TWO IMPULSES 

8.1 Orbital equations in the pure tidal case 

Like the first impulse, to derive the orbit after a second 
impulse we add the component velocity changes to the cur- 
rent velocity, assume the position is unchanged during the 
impulse, and use these quantities to determine a new (rotat- 
ing) p-ellipse orbit. In Section 2 we were able to use the fact 
that the form of the p-ellipse velocity ratio is generally the 
same as that of a circular velocity plus the tidal impulse, to 
derive a very simple orbital solution. With little more effort 
this can also be done after a second impulse. The velocity 
equation ratio analogous to equation (10) can be written, 

(| + S) men sin (mOu) _ v r _ v r3 
[1 + en cos (men)] v g 

Vr2 + AjjV^ COS (2 (02 ~ </>c)) ,„„, 

^ 2 [1-A^sin(2(0 2 -0 C ))]- 1 ' 

The terms on either side of the first equality represent the 
p-ellipse orbit immediately after the second impulse. It has 
an eccentricity en, and its angular coordinate is 9n(t). The 
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final equality gives the immediate post-impulse radial ve- 
locity in the numerator, and the azimuthal velocity in the 
denominator. 

Henceforth, we use the subscripts 0, 1, 2, and 3 to de- 
note the radius and azimuthal positions and the correspond- 
ing velocity components at the following times: just before 
the first impulse, immediately after the first impulse, just 
before the second impulse, and immediately after the second 
impulse, respectively. We use subscript I for parameters of 
the first p-ellipse orbit (between the two impulses), and sub- 
script II for p-ellipse parameters after the second impulse. 
The angle <f) c gives the direction of the second impulse, rel- 
ative to the direction of the first impulse, which is at zero 
azimuth. The amplitude A'jj will generally have a form like 
that given by equations (7) and (9), but with a different 
value. 

Because the orbits are not circular before the second 
impulse, the numerator of the last term in equation (38) 
has a term, v r 2, which is absent in equation (10). However, 
obtaining a solution like that of equations (11) - (13), is not 
much more complicated. We can begin with a very similar 
form for the initial value of the p-ellipse azimuth after the 
second impulse, 0ii o , 



mOuo = 2 (0 2o - 4>c) + ^ + i>- 



(39) 



Then we can use the conditions that the numerators and 
denominators (radial and azimuthal velocities) in equation 
(38) are equal to get, 



A' 



en 



{\+8)m 
'1 



+ 8) m 2 sin 2 (2(^ 2 -0c)) + 



cos(2(0 2 -0 c )) + -^- 



, (40) 



and, 



tan (2 (02 - 4>c) +ip) = 



(§+<5)msin(2 (0 2 - e )) 
[cos (2 (02-^)) + ^]' 



(41) 



like equations (12) and (13). Note that the extra terms from 
equation (16) could also be included in equation (40) with 
the expectation of similar gains in accuracy. 

Like all overlapping wave systems, we expect the or- 
bital phases produced by the two impulses, to generate an 
interference pattern. This complicates the structure of caus- 
tics, for example, or at least generates perturbations to the 
stronger ones. Nonetheless, we will see in the next few sub- 
sections that we can derive some interesting results from the 
solution above. 



8.2 Amplitude of the second impulse 

We continue by looking at the amplitude A' n in more detail. 
This amplitude is an initial condition of the second p-ellipse. 
Following the example of equations (7) and (9) we can write, 



A'n 



12 



2GM C 



12V2V4,, 



_ f 92 



(42) 



Again we use the subscript 2 to denote values at time i 2 
just before the second impulse. For example, V2 is then the 
relative velocity, r m i n 2 is the radius of the second close ap- 
proach, and 52 (= r 2 ) is the radius of the representative star 
at that moment. 

Define the factor 



k\ = 



Ai 



(43) 



which compares the result of changes in the collision param- 
eters r m i n , V , or even M c (in the case of two companions) 
between the two impulses. If the second impulse is simply an 
approximation to an interaction of finite duration, then we 
expect that ItA ~ 1. If it is due to a second encounter, then 
since the relative velocity and distance of closest approach 
probably decrease due to dynamical friction, ItA > 1. For 
encounters with two different companions, values of greater 
or less than one are possible. With the k,A factor equation 
(43) can be used to eliminate An in favor of Ai in Equation 
(42). 

None of the factors in equation (42) depend on azimuth 
except §2, the initial radius at the time of the second im- 
pulse. This dependence is given by the p-ellipse equation 
(equation (1) or (2)). I.e., 



92 



92 qi 



Vmin.2 <?1 ?*min2 



1 + ei cos ( ra0ji ) 



gi [ l-^sin(20 o ) 

Tminl [ 1 + A'j COS (mdl2 ) _ 



(44) 



Equation (10) was used to obtain the last equality. We note 
that 812, for example, is defined as the azimuth of the first 
p-ellipse at the time just before the second impulse. 

Substituting equations (44) into equation (42) and sim- 
plifying yields, 



A'n = k A A! 



1 - J 4jsin(20 o ) 
1 + A'j cos (m6i2 



-,a+s) 



= k' A A'j 



A', sin(20 o ) 



qi \ I r min i 

(*+«) 



1 +A'j cos (m.612) 



(45) 



where the last expression uses equation (9) and the defini- 
tion, 



k' A = k A 



(46) 



The maximum (minimum) value of A' n is obtained 
when sin(20 o ) = — 1(1), cos(m6 l /2) = —1(1), and ej « A' u 
(equation (12)). Those extremal values are, 
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A-Hmax — k A Aj 

and, A IImin — k A Aj 



1 + A'j 
I -A'j 

1 -A', 
1 + A'j 



(HO 



(47) 



In the case where k' A « 1, and is small (e.g., < 0.1), 
then A' IImax « Aj. 

Aj 7 is the net amplitude after the two impulses, so this 
result and the corresponding result for n low amplitude en- 
counters, A' n max « ^'1(1 + 2nA[) confirm the conclusion 
that it would take a very prolonged encounter with a low 
mass companion to significantly perturb the primary. 

When A'j exceeds values of about 0.2 then the values 
of A' IImax begin to grow more nonlinearly. For values of A'j 
approaching 0.5, A' IImax approaches or exceeds unity de- 
pending on the value of 5. If the approximation of equation 
(16) was used to obtain A'j, then the amplification would be 
even greater. Thus, the second impulse can drive a very non- 
linear increase (or decrease) in the perturbation amplitude, 
for orbits with the appropriate phases. 



radius of a stellar orbit can grow geometrically with succes- 
sive impulses, or prolonged disturbances. The dependences 
on S show the dependence on the form of the potential. I.e., 
the effect is roughly doubled in Keplerian potentials relative 
to flat-rotation-curve potentials. 



8.4 Nonlinear tidal disturbances 

In this subsection we extend the mathematical results of Sec- 
tion 5.1 for nonlinear impulses to the case of two impacts. 
As in the tidal case, consideration of multiple impacts intro- 
duces additional parameters, e.g., the amplitude and impact 
parameter of the second interaction. We will not extend the 
sample disc evolution investigations of Section 5.2, since ex- 
ploration of this larger parameter space is beyond the scope 
of this paper. The goal here is primarily to record the equa- 
tions for future use, either with the disturbance functions of 
equations (31) and (32), or with other forms. 

We begin with the velocity matching equations, analo- 
gous to equation (48), with the equations for the radial and 
azimuthal velocities, but written separately rather than in 
a ratio, 



8.3 Maximal excursions 

Equally interesting is the largest radial excursion of the stel- 
lar orbit following a second impulse. We begin by noting that 
equations (17) and (18) apply to the second p-ellipse (gen- 
erated by the second impulse), as well as the first. In par- 
ticular, in the notation of this section, the apoapse radius is 
given by, 



' ap 
<12 



1 + A'jj 

1-A' 



(48) 



with q2 and A' n given by equations (44) and (45), and as- 
suming e/j « A'j j. (The more accurate approximation of 
equation (18) would clutter the already complex equations 
below, without qualitatively changing the conclusions.) 

Under the resonant conditions that led to equation (47), 
and when k' A = 1, equation (45) yields, 



A'n = A'j 



1 + A'j 

1-A'j 



(HO 



(49) 



and 



(1 - A'j)^ +S ) + A'j (1 - A'j)^ +S ) 



(*+«) 



_{1-A'j)^ +5 ) - A'j (1 - A'j)^ +S ) _ 

' l + A^. \(g+«) 

.I" A'j) 

To first order in the amplitude, the apoapse radius is, 



' ap 



1 + 4(1+^, 



(50) 



(51) 



which is about twice the small amplitude increase given by 
equation (18) (at the same level of approximation). 

This shows that in optimal circumstances, the apoapse 



1 

2 

Vr2 



+ S J men sin (mOz) = 

+ A'jj (COS (2 (02 - 0c)) - 7)11 COS (02 - 0c)) , (52) 



ejj cos (77163) — —2rju cos (02 - 

-A'jj (Sin (2 (02 -0c)) 



0c) + VII 

?7//sin(02 - C )) • (53) 



where the subscript notation follows the conventions of Sec- 
tion 4. That is, subscript 2 denotes the values of quantities 
at a time just before the second impulse, subscript 3 denotes 
quantity values immediately after the second impulse, and 
subscript c is used for the azimuth of the companion galaxy 
at the time of the second impulse (or closest approach) . The 
subscript 77 denotes parameters of the new p-ellipse orbit 
after the second impulse, e.g., the eccentricity. The quantity 
rju equals the radius of the star at the time of the second 
impulse divided by the minimum approach radius of the 
companion. As in Section 5 this is the nonlinearity parame- 
ter. 

We assume that we know the parameters of the first 
p-ellipse orbit, generated by the first impulse, and the time 
between impulses (and also the azimuth C ). Then we know 
the subscript 2 quantities, i.e., the positions and velocities 
to which the stars have advanced along their p-ellipse or- 
bits during the time between impulses. In particular, the 
velocity ratio in equations (52) and (53) can be written in 
terms of the stars azimuth before the second impulse and 
the eccentricity of its (first) p-ellipse orbit (as in equation 
(4)), 



v r 2 (| + S) me j sin (17182 
V02 1 + ej cos (mQi) 



(54) 



As in Section 4 we develop the solution for the second p- 
ellipse orbit by adopting the following expression for its ini- 
tial azimuth (see equation (39)), 
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m0 3 = 2 (0 2 - <t>c) + 2 + i>- (55) 

Then, the eccentricity of the new p-ellipse orbit can be de- 
rived by squaring and adding the two velocity equations (52) 
and (53) to eliminate the variable # 3 . The result is, 



V 2 J m e " = 

-1 2 

— + A' n (COS (2 (02 - 0c)) - J?I/ COS (02 - C )) 
_«02 

+ [ - 2//// COS (02 - 0c) + 

- A' u (sin (2 (02 - C )) - r/// sin (0 2 - C )) f ■ (56) 

An equation for the phase -0 can then be derived from equa- 
tion (52) after substituting equation (55), 



COS (2 (02 - 0c) + 1p) = 

+ A' n (cos (2 (02 - C )) - VII COS (02 - 0c)) 
(I + 5) me// 



(57) 



and where equation (54) can be used for the velocity ratio. 
These equations complete the second impulse solution. 

Two impulse models may not accurately represent in- 
teractions of significant duration, but they do illustrate some 
of the qualitative forms that can be generated in such en- 
counters, as discussed in Section 4. The equations of this ap- 
pendix can be readily generalized to describe a sequence of 
nonlinear encounters, to create more realistic semi-analytic 
models of prolonged encounters. Alternate forms can also be 
used in the right hand sides of equations (52) and (53) to 
represent disturbances different than the impulsive Plummer 
potential types used here. 
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